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We consider a tapping dynamics, analogous to that in experiments on granular media, on spin 
glasses and ferromagnets on random thin graphs. Between taps, zero temperature single spin flip 
dynamics takes the system to a metastable state. Tapping, corresponds to flipping simultaneously 
any spin with probability p. This dynamics leads to a stationary regime with a steady state energy 
E{p). We analytically solve this dynamics for the one dimensional ferromagnet and ±J spin glass. 
Numerical simulations for spin glasses and ferromagnets of higher connectivity are carried out, in 
particular we find a novel first order transition for the ferromagnetic systems. 



PACS numbers: 05.20, 75.10 Nr, 81.05 Rm. 

Complex systems such as granular media and spin 
glasses have an exponentially large number Nms of 
metastable states, also called blocked or jammed con- 
figurations, i.e. Nms = exp(A^s_Edio), where SEdw is the 
Edwards entropy per particle ||l[ . In granular media the 
thermal energy available is not sufficient to allow the re- 
arrangement of a single particle and hence the system 
is effectively at zero temperature in the thermal sense 
and when not perturbed are stuck in a metastable state. 
Edwards has proposed 0] that if one slightly perturbs 
such systems by external forcing such as tapping, the 
asymptotic measure over the metastable states satisfy- 
ing the appropriate macroscopic constraints is flat, this 
idea has recently been further investigated In a, 

now classic, experiment on granular media B dry hard 
soda glass spheres placed in a glass tube are tapped by 
using a piston to move the tube vertically through a sine 
cycle. The tapping parameter F is defined to be the ra- 
tio of the maximal acceleration due to the piston in the 
cycle to g the acceleration due to gravity. After an ini- 
tial irreversible curve obtained by increasing the tapping 
amplitude slowly, the system arrives on a reversible curve 
where the density is a monotonic function of F, the high- 
est packing densities being obtained at lowest tapping 
rate. Numerical simulations on granular media reveal 
similar behavior. There has been a considerable theoreti- 
cal effort to understand the dynamics of granular media, 
especially the slow relaxation towards the final asymp- 
totic density under tapping, both theoretically and by 
simulations on various models |^-^. Our aim is to gain 
some understanding of the stationary regime of tapped 
systems such as granular media. 

The possibility of using spin glasses as a paradigm for 
granular material was first suggested in , and recently 
it has been shown in ||l^,|ll| that spin glasses and ferro- 
magnets on random thin graphs have an extensive Ed- 
wards entropy with respect to single spin flip dynamics. 
In this paper we study a very natural tapping dynamics 
on these systems: We allow the system to evolve un- 



der a random zero temperature single spin flip dynamics 
where only moves which reduce the energy are allowed 
(falling). When the system is blocked we tap it with 
strength p £ [0, 1/2], that is to say each spin is flipped 
with a probability p, the updating at this point being 
parallel. The system is then evolved by the zero tem- 
perature dynamics until it becomes once again stuck, the 
tapping is then repeated. Physically this corresponds to 
assuming that in granular media the relaxation time to 
a new metastable state is much shorter than the time 
between taps. The same tapping dynamics has also been 
introduced independently in the context of three spin 
ferromagnetic interactions on thin hypergraphs also 
in the goal of studying the dynamics of granular media. 
We find that a stationary regime is reached after a suf- 
ficiently large number of taps, characterized by a steady 
state energy E{p) (analogous to the stationary density 
- the same analogy as used in [|j). We then develop 
a mean field theory for the dynamics under falling then 
tapping, which appears to be exact in the case of the 
above mentioned one dimensional system and one may 
calculate E{p) within this scheme, the results being in 
perfect agreement with the numerical simulations. We 
also examine the tapping of spin glasses and ferromag- 
nets of higher connectivity. In the case of the ferromag- 
net we find numerically that there exists a critical value 
Pc of p such that for p > pc, E(p) > Ecs where Egs 
is the energy of the ground state and the inequality is 
strict, and that for p < pc , E{p) — Eqs^ hence in the 
ferromagnetic system there is a first order phase transi- 
tion with the tapping dynamics (in contrast to the usual 
thermodynamic ferromagnetic transition in these systems 
which is second order |l^]). The spin glass is clearly far 
from a realistic realization of a granular media, however 
the fact that it has extensive entropy of blocked states 
and the obviously natural form of the tapping dynamics 
implemented makes it a natural testing ground for ideas 
about dynamics and possible thermodynamics of systems 
such as granular media. 
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A random thin graph of connectivity c is a collection 
of N points, each point being linked to c of its neighbors. 
The spin glass/ferromagnet model we shall consider has 
the Hamiltonian 



H = - 



1 



2 'Jij'Hj^i'-'j 



(1) 



where the Si are Ising spins, is equal to one if the 
sites i and j are connected. In the spin glass case the 
are taken from a binary distribution where Jy — —1 with 
probability half and = 1 with probability half. In the 
ferromagnetic case Jy = 1. The definition of the total 
number of metastable states in this system is the number 
of states where any single spin flip does not increase the 
energy of the system. 

The one dimensional spin glass/ferromagnet: We re- 
mark that by a gauge transformation the one dimensional 
ferromagnet and ± J spin glass are equivalent and place 
ourselves, for transparency, in the context of the ferro- 
magnet. In the initial configuration, we take the prob- 
ability that a given spin is different to its left neighbor 
to be a. Hence if a = we have an initially ferromag- 
netic configuration, if a = 1 it is an antifcrromagnctic 
configuration, the case a — 1/2 corresponds to a com- 
pletely random configuration. The initial energy per spin 
is Eo = -1 + 2a. 

For a given site define x to be the difference between 
the number of unsatisfied and satisfied bonds. Hence x is 
the local field on the spin; if a; > then the spin can flip 
bringing about the change x — > —x. Denote by P{x, k) 
the probability that the site of interest has local field 
x after a total of k attempted random single spin flips 
during the falling dynamics. We define /+ and /_ the 
probabilities that a neighboring spin can flip conditional 
on the fact that the bond with the site we are interested 
in is not satisfied or satisfied respectively. Given that a 
site has local field x, it must have (c -I- x)/2 unsatisfied 
bonds and {c—x)/2 satisfied bonds; using Bayes' theorem 
we therefore obtain 



f± = 



j:^pix)ic±x) 



(2) 



For the spin considered the possibilities between time 
k and fc -I- 1 are 

• X > 0. Then the spin can flip and x goes to —x; 

• A neighboring spin has a positive local field and so 
can flip. In that case, x goes to a; + 2 or to x — 2; 

• Neither the spin considered nor any of its neighbors 
flips, and so x stays x. 

Assuming that the distribution at every site is given by 
P{x, k) and assuming independence between the values 
of X from site to site (the mean field approximation) we 



obtain, taking the limit N oo and introducing the 
continuous time t = k/N, 



dP{x) 
dr 



= e{-x)p{-x) + e{-~x)p{x) - (c + i)P{x) (3) 

+ P{x) + 



The case where c = 2 is accessible to analytic solu- 
tion. Defining: u(t) = P(-2,r), v{t) = P(0,t) and 
'w{t) = P(2,t), the solution to Eq. (pi) is given by 



w{t) — w{0) exp — r + 



A(0) + 2 

v{t) = -2u;(t) + w{0){X{0) + 2) exp 



(4) 



A(0) + 2 



(e--l) 



where A(0) = v{0)/uj{0). Consequently after the sys- 
tem has fallen into a metastable state we find the final 
values v{oo) = (u(0) + 2w{0)) e~^^s^ , uj{oo) = and 
m(oo) = 1 — v(oo). 

For initial conditions used here: u(0) = (1 — a)^, 
v{0) = 2a(l — a) and w{0) = o/ . Hence i?/, the aver- 
age energy of the metastable state into which the system 
falls, is given by 



= + vice) = -1 + 2ae" 



(5) 



This result can be shown to be exact by a combinatorial 
calculation and has also been checked in our simulations 
p3| . For the completely random initial configuration, 
where a — 1/2, Ef = —0.632121 is in fact maximal 
This calculation demonstrates two important points: 

• The final value of the energy Ef depends strongly 
on the initial configuration, in addition E^ is not a 
monotonic function of Eq. This means that con- 
figurations of higher initial energy can fall into 
metastable states of lower energy than initial con- 
figurations of a lower energy. 

• The system does not fall into a state of energy cor- 
responding to the maximum of NMsiE): the to- 
tal number of metastable states of energy E per 
spin. In (l^Jll[] it was shown that Nms{E) ^ 
exp{NsEdw{E)) where SEdw{E) is a concave func- 
tion peaked at E* ^ -1/V5 « 0.44721. Hence 
even if the total number of metastable states is 
dominated (in the thermodynamics sense) by those 
of energy E* , generic initial conditions always seem 
to lead to an energy lower than this |2||~ 

Tapping the system with probability p, starting from 
the values {u(cxd), w(oo), w(cxd)}, we obtain the new 
tapped values {u' {0) , v' {0) , w' {0)} . Defining q = {1 —p), 
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the relations between the old and tapped probabilities 
are: 



u'(0) = (1 — 3pq) u(oo) + pq v{co) 
v'{0) = 2pq u{oo) + (1 - 2pq) v{oo) 
w'{0)^pq 



(6) 



After another zero temperature evolution of the system, 
it reaches a new local energy probability distribution with 
w'{oo) = 0, v'{oo) and m'(oo) = 1 — v'{co). The fixed 
point of this dynamics is given by 



Vs{p) = i-ipq + Vs{p){l- 4pq)) exp 



Apq 



Apq 



Vs{p){l - 'ipq) 
(7) 



the subscript s indicating steady state. We remark that 
this solution contains one of the main features of our 
numerical simulations, that is reversibility. Here the 
asymptotic distribution is independent of the initial con- 
ditions and depends only on p. Eq. can be solved 
numerically and the result is shown in Fig. (1) in com- 
parison with the numerical simulations which we see is 
excellent. The small p behavior of E{p) from is: 
E{p) ^ -l + ^+0{p). 

Systems with c > 2: cases. In [|lO|,0 an an- 
nealed approximation to the Edwards entropy per spin 
of metastable states at fixed energy E, SEdw{E) — 
ln{{NMs{E))) /N was carried out for spin glasses and 
ferromagnets on random thin graphs. There is an energy 
threshold E* above which the results are the same for 
the ±J spin glass and the ferromagnet; below E* the 
ferromagnet has more metastable states (which aquire a 
non-zero magnetization pl]|). Hence, as far as metastable 
states are concerned, both ferromagnet and spin glass are 
the same above E* , that is the effect of loop frustration 
is negligible. In this regime, one suspects that the zero 
temperature dynamics are the same. In particular, nu- 
merical simulations with 100 samples of = 10000 sites 
for connectivities of 3, 4 and 5 have found the same Ef 
for the spin glass and ferromagnet to very good accuracy 
(the relative error is about 10~^). The results of tap- 
ping experiments on systems with c = 3 are displayed in 
Fig. (2). There is a critical tapping rate Pc above which 
the curves of E{p) versus p are the same for the spin 
glass and ferromagnet. At pc ~ 0.249 the ferromagnet 
undergoes a phase transition such that for p < pc, the 
steady state reached is the ground state. Finite size ef- 
fects have been studied and reveal that the transition 
appears to be first order (in as far that E{p^) ^ E{p'^)). 
As shown in Fig. (3), there is a coexistence of two phases 
at Pc] one sees two separated peaks in the distribution of 
the internal energy and not a single peak as one would 
expect for a second order transition. In [ pT| it was shown 
that for the ferromagnet the annealed approximation to 
the Edwards entropy as a function of E is concave for 



E > E* and convex for E < E*. The value of £'(p+) 
estimated from the tapping experiments are very close 
to those obtained for E* in the energy at which 
SEdw{E) becomes convex (for c = 3, i?* = —1.0714 an- 
alytically and E(j)+) « -1.075 ± 0.005 from the sim- 
ulation). Encouraged by this observation we will try to 
make a tentative link with a possible thermodynamics for 
such systems. We consider a partition function inspired 
by the flat Edwards measure over metastable states [^] 
Z = / dENMs{E) exp{—f3NE) where /3 is a Lagrange 
multiplier corresponding to the inverse Edwards temper- 
ature which depends solely on p, and not on i?, and is a 
monotonically decreasing function of p for p S [0+, 1/2]. 
The monotonicity hypothesis is supported by the sim- 
ulation results that E{p) decreases with decreasing p. 
Clearly the energy which dominates in the sum is that 
obeying ^^"^g]^^^^ — P — 0. However if this saddle point 
gives a true maximum of the action the Edwards entropy 
must be concave for the energy considered to be thermo- 
dynamically stable. Hence one is lead to conclude that 
for E < E* the only stable energy is the ground state. 
Finally we remark for c > 3 that the mean field equa- 
tions of the previous section may be numerically solved 
and give reasonable results for E{p) for p > Pc |0 , which 
again reinforces the assertion that for higher energies the 
spin glass and ferromagnetic dynamics are the same (the 
mean field theory does not distinguish between the two).. 

No observable singularities are seen in our simulations 
of spin glasses, but again we found that for all the sys- 
tems studied that E{p) decreased as a function of p. As 
p ^ we found that generically E{p) ^ £^(0+) + Ap^ . 
For the one dimensional case we know analytically that 
6 = 1/2, however for c = 3, 4, and 5 we found that 
6=1. We also carried out simulations on the totally 
connected Sherrington Kirkpatrick ± J spin glass for sys- 
tems of size 400, 600 and 800 spins where an average was 
taken over 1000 samples. Here we used a deterministic 
sequential (rather than random) update for the zero tem- 
perature dynamical evolution - as used by Parisi in 
where the numerical evaluation of Ef was undertaken. 
Interestingly we found that the finite size scaling intro- 
duced by Parisi for Ef worked extremely well for E(p) 
that is E{p, oo) ^ E{N,p) - 1/4A^3, for all p. In ad- 
dition we found that for p small E{p) — i?(0+) ~ Ap^ 
with 6* « 0.4 ± 0.1 |l3|. We emphasize that all of the 
numerical simulations found here (for spin glasses and 
ferromagnets) exhibit perfect reversibility for E{p). In 
addition the values oi Ef found in all systems were less 
than the value of the energy dominating thermodynami- 
cally in SEdw{E) (calculated in JlO| , pl| ), showing that the 
falling dynamics in strongly non erg odic §. 

Granular media are examples of systems having an ex- 
tensive entropy of metastable states. In such systems the 
role of thermal fluctuations are negligible and in order to 
evolve one must apply some external tapping mechanism. 
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One would ultimately like to be able to formulate some 
sort of thermodynamics for such systems. The propo- 
sition of Edwards j|] for such a thermodynamics is an 
important step in this direction and has had some suc- 
cess but in the same paper it was shown not gener- 
ically to be true. A more general understanding of the 
asymptotic states of tapped systems has also far reaching 
implications for computer science as the tapping mech- 
anism studied here is similar to certain algorithms used 
in optimization problems. We have presented what ap- 
pears to be the exact solution to the problem of tapping 
a one dimensional spin glass or ferromagnet. The fixed 
point equation for Eq. (|^) may have a thermodynamic in- 
terpretation which presents an open challenge to find. In 
a wide context of models we confirm the observations of 
iQ-^, that if one reduces the strength of tapping, then the 
compaction process, corresponding here to the reduction 
of the energy of the system, becomes more efficient. The 
existence of a first order type phase transition for tapped 
ferromagnets on random thin graphs is of great interest, 
the possible explanation using the calculations of |0 on 
the Edwards entropy for this system indicates the possi- 
bility of constructing a general theory for the thermody- 
namics and phase transitions in tapped systems. 
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Figure Captions 

Fig. 1. Comparison between numerical simulations of tapping experiments (b) and the analytical result (a) obtained 
with (0). 

Fig. 2. Numerical simulations of tapping experiments for the spin glass (c) and the ferromagnet ((a), (b) and (d)) for 
c = 3 for iV = 1000 ((c) and (d)), N = 2000 (b) and N = 10000 (a). The inset shows the scaling E^^^ = f{N{p-p^)) 
for p~pc for N = 400, N = 1000 and N = 2000. 

Fig. 3. Histogram of the distribution p{E) of energy during a tapping simulation at p w for a single run of 50000 
taps and N — 5000 spins. 
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